Sparling two-forms, the conformal factor and 
the the gravitational energy density of the 
teleparallel equivalent of general relativity 

J. W. Maluf* 
Departamento de Fisica 
Universidade de Brasilia 

CP. 04385 
70.919-970 Brasilia, DF 
Brazil 

Abstract 

It has been shown recently that within the framework of the telepar- 
allel equivalent of general relativity (TEGR) it is possible to define in 
an unequivocal way the energy density of the gravitational field. The 
TEGR amounts to an alternative formulation of Einstein's general rel- 
ativity, not to an alternative gravity theory. The localizability of the 
gravitational energy has been investigated in a number of spacetimes 
with distinct topologies, and the outcome of these analises agree with 
previously known results regarding the exact expression of the grav- 
itational energy, and/or with the specific properties of the spacetime 
manifold. In this article we establish a relationship between the ex- 
pression for the gravitational energy density of the TEGR and the 
Sparling two-forms, which are known to be closely connected with 
the gravitational energy. We also show that our expression of energy 
yields the correct value of gravitational mass contained in the confor- 
mal factor of the metric field. 
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I. Introduction 

The problem of a consistent definition of the gravitational energy has 
been addressed over the years by means of rather different methods. The 
traditional approach amounts to defining the total gravitational energy of 
an asymptotically flat spacetime through the use of pseudotensors. A closely 
related approach consists in associating the gravitational energy with the sur- 
face term that appears in the total Hamiltonian of the gravitational field. An 
alternative way of defining the gravitational energy is provided by the con- 
cept of quasi-local energy. The quasi-local definition of energy, momentum 
and angular momentum associates these quantities with an arbitrary two- 
surface S of an arbitrary spacetime manifold. This approach is presented in 
ref . |]TJ , where a comprehensive list of references on the problem of energy in 
general relativity is presented. 

A recent approach to the problem of the definition of the gravitational 
energy has arisen in the framework of the teleparallel equivalent of general 
relativity (TEGR)JJ, [3], |J. The latter is an alternative formulation of general 
relativity. This formulation is established by means of the tetrad field e afl and 
the spin connection tu^ a b, which are totally independent field quantities, even 
at the level of field equations. The metric tensor obeys Einstein's equations. 
However, the action integral is constructed in terms of the torsion tensor. 
The Lagrangian density of the TEGR is a functional of the torsion tensor, 
but is precisely equivalent to the ordinary scalar density eR(e aiM ), provided 
the curvature tensor Rab^i^) vanishes. In fact the vanishing of R a b^u{^) is 
one basic motivation for considering the TEGR, since this property leads to 
the establishment of a reference space||. 

The definition of the gravitational energy density emerges in the context 
of the Hamiltonian formulation of the TEGR. The latter is considered in 
ref.||. It has been shown that, under a suitable gauge fixing, the Hamiltonian 
formulation of the TEGR is well defined, as the constraints turn out to be 
first class. The major property of the Hamiltonian formulation is that the 
Hamiltonian constraint C = can be written as 

C = H — Eadm = , 

in the case of asymptotically flat spacetimesQ, [|. H is interpreted as the 
effective Hamiltonian. The Arnowitt-Deser-Misner (ADM) energy|| (Eadm) 
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arises in the Hamiltonian constraint upon integration over the whole three- 
dimensional spacelike hypersurface of a scalar density e(x), which can be 
written as a total divergence of the form g^<9j(eT J ), where T % is the trace 
of the torsion tensor. We propose that e(x) represents the gravitational 
energy density for spacetimes with any topology. In fact we have applied it 
to the determination of the distribution of gravitational energy in de Sitter 
space and found that the resulting analysis is in total agreement with the 
phenomenological features of this space 0. Therefore we assert that the 
Hamiltonian constraint equation can be written as C = H — E = for any 
spacetime. 

The existence of the gravitational energy density as given above allows 
the conclusion that the gravitational energy is localizable, in the same way as 
the electromagnetic energy is also localizable. The very idea of a black hole 
lends support to this conclusion. If there exists an ammount of gravitational 
energy (mass) inside the event horizon of a black hole, then by means of a 
coordinate transformation we do not expect to make this energy vanish. 

In this article we establish a relationship between the gravitational energy 
density e(x) and the Sparling two-forms. The latter are known to yield the 
gravitational energy and momentum of an asymptotically flat spacetime by 
means of an integration over a two-surface at infinity. In fact the Sparling 
forms are closely connected with various pseudotensors of energy-momentum 
of the gravitational field. We also show that under integration e(x) captures 
the information about the gravitational mass contained in the conformal 
factor of the metric field, provided the appropriate boundary conditions at 
infinity are imposed. Thus our expression of energy encompasses the features 
of distinct approaches to the definition of the gravitational energy. 

We remark that although there exist in the literature several expressions 
for the definition of the total gravitational energy, the present expression 
is unique as it is the only one that possesses the feature of localizability. 
It is precisely this feature that allows obtaining the striking result regarding 
rotating black holes (cf. ref. 0), namely, that the energy contained within the 
outer horizon of a Kerr black hole is essentially equal to twice its irreducible 
mass. 

Notation: spacetime indices fx, p, ... and local Lorentz indices a, b, ... run 
from to 3. In the 3+1 decomposition latin indices from the middle of 
the alphabet indicate space indices according to /x = 0,i, a = (0),(z). 
The tetrad field e a M and the spin connection u^b yield the usual definitions 
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of the torsion and curvature tensors: R a b^u = d^Uy a j, + a c uo v c b — 
T a = d^e a u + u> fl a be b u — • ••• The flat spacetime metric is fixed by 
V(o)(o) = -I- 



II. The TEGR 

Throughout the paper we will consider only asymptotically flat space- 
times, because it is in this geometrical framework that the Sparling two-forms 
are relevant for energy considerations. The Lagrangian density of the TEGR 
in empty spacetime is given by 

L(e,u,X) = -ke{ l -T abc T abc + l -T abc T bac - T a T a ) + e\ ab ^ R ab , u {u) . (1) 

where k = G is the gravitational constant; e = det(e a M ), \ abtil/ are 

Lagrange multipliers and T a is the trace of the torsion tensor defined by 
T a = T b ba- The tetrad field e a ^ and the spin connection u>^ a b are completely 
independent field variables. The latter is enforced to satisfy the condition of 
zero curvature. Therefore this Lagrangian formulation is in no way similar 
to the usual Palatini formulation, in which the spin connection is related to 
the tetrad field via field equations. 

The equivalence of the TEGR with Einstein's general relativity is based 
on the identity 

eR(e,u) = eR(e) + e{-J abc T abc + ^T abc T bac - T a T a ) - 20 Ai (eT") , (2) 

which is obtained by just substituting the arbitrary spin connection u^ a b = 
°^nab{e) + K^ a b in the scalar curvature tensor R(e,u) in the left hand side; 
°^V&(e) is the Levi-Civita connection and K^ ab = \e a x e b v iT\ ViV + T uXfl - 
T,j, v \) is the contorsion tensor. The vanishing of R a which is one of the 

field equations derived from (1), implies the equivalence between the scalar 
curvature R(e), constructed out of e a M only, and the quadratic combination 
of the torsion tensor. It also ensures that the field equations arising from 
the variation of L with respect to e a M are strictly equivalent to Einstein's 
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equations in tetrad form. Let = denote the field equations satisfied by 
e aA \ It can be shown by explicit calculations that 

^ = \{Ra,{e) - \e a ,R{e)} . (3) 

We refer the reader to refs. |3], |3J] for additional details. 

In the Lagrangian density (1) we have not included the total divergence 
d^(eT^). The reason is that the variation of the latter in the action integral 
causes the appearance, via integration by parts, of surface terms that do not 
vanish when r — > oo. In this limit we should consider variations in g^ u or 
in e afJi that preserve the asymptotic structure of the flat spacetime metric; 
the allowed coordinate transformations must be of the Poincare type. The 
variation of <9 M (eT M ) at infinity under such variations of e afJ- does not vanish. 
Therefore the total divergence has to be dropped down. On the other hand, 
all surface terms arising from partial integration in the variation of the action 
integral constructed out of (1) vanish in the limit r —>■ oo. The action based 
solely on the quadratic torsion terms does not require additional surface 
terms, as it is invariant under transformations that preserve the asymptotic 
structure of the field quantities. This property fixes the action integral, 
together with the requirement that the variation of the latter must yield 
Einstein's equations, which is actually the case in view of (3) (the Hilbert- 
Einstein Lagrangian requires the addition of a surface term for the variation 
of the action to be well defined; a clear discussion of this point is given in 
ref.H). Unfortunately in refs. [2-5] we have not given the correct argument for 
the suppression of d^eT^). We argued, erroneously, that the variation of this 
divergence vanishes at infinity when we consider variations of e afl as defined 
above. This divergence has to be subtracted by hand, in the same way that 
the Hilbert-Einstein Lagrangian density requires an additional surface term. 

In what follows we will be interested in asymptoticaly flat spacetimes. 
The Hamiltonian formulation of the TEGR can be successfully implemented 
if we fix the gauge uo a b = from the outset, since in this case the constraints 
(to be shown below) constitute a first class set||. The condition u 0a b = 
is achieved by breaking the local Lorentz symmetry of (1). We still make 
use of the residual time independent gauge symmetry to fix the usual time 
gauge condition e^) = e( )j = 0. Because of uo a b = 0, H does not depend 
on p kab ^ the momentum canonically conjugated to Ukab- Therefore arbitrary 
variations of L = pq — H with respect to P kab yields Cj kah = 0. Thus in view 
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of uoab = 0, Ukab drops out from our considerations. The above gauge fixing 
can be understood as the fixation of a global reference frame. 

As a consequence of the above gauge fixing the canonical action integral 
obtained from (1) becomes [§J 

A TL = J d 4 x{Tl^ k e m - H} , (4) 



H = NC + N l C t + E mn U mn + J—d k (NeT k ) + d k {W k N,) . (5) 

07TG 

N and N* are the lapse and shift functions, and E mn = — S„ m are Lagrange 
multipliers. The constraints are defined by 

C = d&keT 3 ) - keX kij T kij - -^(H"1I ; , - ^n 2 ) , (6) 

C k = -e U)k d t U^ - U^T (j)ik , (7) 

with e = det{eu-\ k ) and T % = g tk e^ l T^i k . We remark that (4) and (5) are 
invariant under global SO (3) and general coordinate transformations. 

We assume the asymptotic behaviour e^ k ~ rjj k + \hj k {-) for r — > oo. 
In view of the relation 



— - / (PxdjieT*) = — — I dS k (d t h lk -d k h u ) = E ADM (8) 
8nG J IQttG Js 

where the surface integral is evaluated for r — > oo, the integral form of the 
Hamiltonian constraint C = may be rewritten as 

J d 3 xheZ kij T kij + 4^(11%* - in 2 ) j = E ADM . (9) 

The integration is over the whole three dimensional space. Given that dj{eT^) 
is a scalar density, from (7) and (8) we define the gravitational energy density 
enclosed by a volume V of the space as 

E -Sfl/v As * Tf) ' (10) 
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It must be noted that E g depends only on the triads e^)i restricted to a 
three-dimensional spacelike hypersurface; the inverse quantities e^ k > 1 can be 
written in terms of e^v. From the identity (3) we observe that the dynamics 
of the triads does not depend on u^. Therefore E g given above does not 
depend on the fixation of any gauge for u^. 

The reference space which defines the zero of energy has been discussed 
in ref.[[|. Here we will briefly present the main ideas about its fixation. 
The establishment of a reference space requires the concepts of holonomic 
and anholonomic transformations of coordinates. Let us consider the Eu- 
clidean space with metric = (+ + +)> which is the spatial section of the 
Minkowski metric. We introduce coordinates such that the line element 
of the Euclidean space is written as ds 2 = rj^^dq^ dq^\ We consider next 
a coordinate transformations dq^' = j(x)dx^ , which allows us to rewrite 
ds 2 = V(i)(j) e m{ x ) e n (x)dx m dx n = g mn dx m dx n . This transformation can 
be holonomic or anholonomic. 

If the relation dq^ = j(x)dx^ can be integrated over the whole three- 
dimensional space, the transformation gW — > x J corresponds to a single- 
valued global transformation, and therefore it is called holonomic. Both sets 
of coordinates, {q^} and {x- 7 }, describe the three-dimensional Euclidean 
space, and we have necessarily that j is a gradient vector, i.e., j = 
TjV- However, in general dq^ = jdx^ cannot be globally integrated, 
since j may not be written as the gradient of a function, namely, j 
may not be of the type e'''j = djq^. If the quantities e^'j are such that 
dj£ k — @k e j 0; then the transformation is called anholonomic. 

For triads which are gradient vectors, the torsion tensor T(j)jfe = dje^k — 
dk&{i)j vanishes identically. A crucial result is that T^jk vanishes if and 
only if {e"j} are gradient vectors ||. In the framework of the Hamiltonian 
formulation of the TEGR the gravitational field corresponds to a configu- 
ration for which T^j k ^ 0. We conclude that every gravitational field is 
anholonomically related to the three-dimensional Euclidean space, which is 
to be taken as the reference space. Since the torsion tensor vanishes for the 
latter, as we have seen, the gravitational energy of the reference space is zero. 



III. The Sparling two-forms and its relation with the TEGR 

The Sparling two- forms a a are defined byfHf 
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o«=-\s M ^Ae\ (11) 



where e a bcd is the totally antisymmetric Levi-Civita tensor such that £(o)(i)(2)(3) 

i; 

by 



1 ; T ab and e a are one-forms, T ab = T M ab dx^ , e a = e a ^dx^, which are related 



de a + r% A e b = 0. (12) 
In components, T^b turns out to be the Levi-Civita connection: 

r^afe — 6 ^(y^abc ^bac ^cab) 

Qabc = e aX (e b u d l/ e c x - e c v d u e b x ) . 

It is known that in a coordinate basis which is asymptotically cartesian the 
Sparling forms are related to various pseudotensorsjll], |12| . In particular, 
it allows a definition of the total energy-momentum P a of the gravitational 



field [|12| 



P a = / v a , (13) 

8n& JdT, 



where <9E actually represents a spacelike surface S at infinity (a> afe of ref. |L2| 
differs by an overall sign from T ab ). We mention that the connection between 
the Sparling forms with the Brown- York expression for quasi-local energy has 
been investigated in ref. |13fl . 

In order to establish the relation between P(p) given by (13) and expression 
(10) we need to rewrite <T(o) in components. We have: 



a (o) = — ^(o)bcdT bc A e d 

= ~e m){k) T^ v dx» /\dx v 
= -^s mm T m ^e^ n dx m Adx n . 

The last equality is obtained in view of the fact that a® is integrated over 
a spacelike surface. 
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Let us introduce the surface element dSi defined by 



dSi = -Eijkdx^ A dx k , 



where e 123 = 1. In view of the relations 



P (0 P (J) P (k) mnl _ (i)(j)(k) 



we find that cr( ) can be rewritten, after a number of manipulations, as 

a (0 ) = -ee« m e^>r nW0 .)^ m . (14) 
Therefore substitution of (14) in (13) leads to 

P ( o) = ^j/ )m e {3)n T nm) dS m = ^j v d m {ee^>T nm) )d*x. 

(15) 

It is not difficult to verify that if we assume the asymptotic behaviour e^j ~ 
Vij + i^y(r) wrien r ~" ^ °°) an d impose the usual time gauge condition 
e(j) = j = (the latter are tetrad components of the four- dimensional 
spacetime), then the expression above yields the ADM energy. 

The relation of (15) with the energy expression (10) can now be estab- 
lished in a straightforward way. The equivalence between the two expressions 
rests on the identity 

d k (ee^ k e^ n r n{m ) = d k {eT k ) , (16) 

which can be verified by just substituting ^ m {i){j) 011 the left hand side of the 
equation above. Note that because of the time gauge condition, which is also 
imposed in the Hamiltonian formulation of the TEGR, r m (i)(j) is constructed 
out of the triads components restricted to the three-dimensional spacelike 
hyper surf ace. 

The equivalence between (10) and (15) can, however, be established only if 
the integration is performed over the whole three-dimensional space. The rea- 
son is that a a defined by (11) is normally considered a non-invariant quantity, 
as it transforms inhomogeneously under e%(x) — * e a ^(x) = A a b(x) e b ^(x), 
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where A a b(x) belongs to the local SO(3, 1). On the other hand, in the frame- 
work of the TEGR <9j(eT') is a scalar density, invariant under global SO(3) 
transformations. This is a necessary requirement in order to arrive at a 
Hamiltonian formulation with only first class constraints. 



IV. The conformal factor of the metric 

An alternative way of defining the energy of an asymptotically flat grav- 
itational field is by identifying it with the 0(r~ l ) part of the conformal 
factor of the metric [II]]. However, this identification is only possible if some 



boundary contitions are imposed. Suppose that the metric field satisfies the 
following conditions in the asymptotic limit r — > oo: 

gij = ^ 4 5ij + hij ; if) = 1 + 0(r _1 ) , 

K = I to* kj = 0(r- 2 ) ; dfiij = 0(r- 3 ) . (17) 

If this asymptotic behaviour is verified, the following definition is suggested 



for the energy of the gravitational field [14 



E = dSidiij. (18) 

Z7TG JS^oo 

We will show that a similar statement can be made in the context of 
expression (10). Let us consider the conformal factor if) as above, with the 
asymptotic behaviour if) = 1 + 0(r _1 ), and write 

e(k)i = ^ 2 °e {k )i + e {k)i , (19) 

where {°e(k)i} are triads of the flat three-dimensional space (in cartesian 
coordinates, °e^)i = Ski) and e^)i is such that for r — > oo we have e^)i = 
0(r _1 ). In terms of these quantities we can construct the trace of the torsion 
tensor T* = g ik e {l)j T (l)jk : 

T = g^ e ^[°e {l)k d 3 ilj 2 - e {l)j d k iP 2 + dfe {l)k -d k ~e m }. 

In order to evaluate (10) and compare it with expression (18) we will have 
to obtain the value of eT l on a surface S at infinity. In this limit, the last 
two terms of T % of the expression above contribute to the surface integral as 
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e^dje^^dji^e^) , (20) 

e^d k e {l)j - d k { °e^e {l)j ) . (21) 

We would like to have these quantities vanishing under integration. For this 
purpose it is necessary to have them falling off as ^ at infinity. This will be 
the case if we require 

dj( e {l)k ) = 0(r~ 3 ) ; °e^e (0 , = 0{r' 2 ) . (22) 

We observe that conditions (22) are the equivalent in triad form of conditions 
(17): the first condition above is equivalent to djhy = 0(r~ 3 ), the second 
one to trg hij = 0{r~ 2 ). Therefore we assume them to hold, together with 
(19), for asymptotically flat gravitational fields. As a consequence we obtain 



1 



dS t eT = -—— i dSM, (23) 



+oo 



8irG Js^oo 2ixG Js~ 

which establishes the equivalence with (18). 

As a simple application of expression (23) let us consider the Schwarzschild 
metric in isotropic coordinates. The metric for the spacelike hypersurface is 
given by 

ds 2 = ( 1 + ^ )\dr 2 + r 2 d6 2 + r 2 sin 2 6 d(f> 2 ) , (24) 



2rJ 

where, for convenience, we make use of spherical coordinates. We have also 
made G — 1. The triads associated with this metric read 

/ 1\) 2 sinOcoscj) r ip 2 cosOcoscf) —r vp 2 sin6sin<p \ 
e(fe)i = ip 2 sin6siri(f) r ip 2 cosOsincf) r vp 2 sin6cos<fi (25) 
\ ip 2 cosO — r ip 2 sinO J 

where ip — (1 + ^). In the expression above (k) and i are line and column 
indexes, respectively. It is easy to verify that ie^)j yields precisely the 
metric components of the line element (24). We will exempt from presenting 
the details of the calculations, which in fact are not complicated. The only 
contribution to the surface integral is given by 
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eT 1 = eg lj e {k)i T { 



(k)ij 



—4 r 2 sinO ip-^- ■ 
or 



We find (recall that G = 1) 




87T /s^oo 



dSeT 1 = — 



d9 d(p sinO 2m 



= m , 



(26) 



as expected. 
V. Comments 

We have seen that expression (10) encompasses the features of previ- 
ous, distinct approaches to the definition of the gravitational energy. Fur- 
thermore, it is still possible to consider localized gravitational energy. In 
this paper we have considered the definitions of gravitational energy (i) con- 
structed out of the Sparling forms, and (ii) constructed out of the confor- 
mal factor of the metric field. None of these approaches arise in a natural 
or conventional way in the Lagrangian or Hamiltonian formulations of gen- 
eral relativity. Moreover, they provide only the total gravitational energy 
of asymptotically flat fields. We have seen that these definitions are in fact 
related to the definition (10) for the gravitational energy, which in turn does 
appear as one element of the Hamiltonian constraint C in the canonical for- 
mulation of the TEGR. Thus expressions (13) and (18) can be understood as 
different manifestations of the energy expression (10), when the integration 
is performed over the whole three-dimensional space. This result supports 
the general validity of expression (10). 
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